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Abstract 

We have performed a detailed investigation of the electrical properties of a graphene monolayer 
which is modulated by a weak one dimensional periodic potential in the presence of a perpendicular 
magnetic field. The periodic modulation broadens the Landau Levels into bands which oscillate 
with B. The electronic conduction in this system can take place through either diffusive scattering 
or coUisional scattering off impurities. Both these contributions to electronic transport are taken 
into account in this work. In addition to the appearance of commensurability oscillations in both the 
collisional and diffusive contributions, we find that Hall resistance also exhibits commensurability 
oscillations. Furthermore, the period and amplitude of these commensurability oscillations in the 
transport parameters and how they are affected by temperature are also discussed in this work. 



1 



I. INTRODUCTION 



Recent successful preparation of a single layer of graphene has generated a lot of inter- 
est in this system as experimental and theoretical studies have shown that the nature of 
quasiparticles in this two-dimensional system is very different from those of conventional 
two-dimensional electron gas (2DEG) systems realized in semiconductor heterostructures. 
Graphene has a honeycomb lattice of carbon atoms. The quasiparticles in graphene have 
a band structure in which electron and hole bands touch at two points in the Brillouin 
zone. At these Dirac points the quasiparticles obey the massless Dirac equation. In other 
words, they behave as massless, chiral Dirac Fermions leading to a linear dispersion relation 
Ek = hvpk (with the characteristic velocity vp — lO^m/s). This difference in the nature of 
the quasiparticles in graphene from a conventional 2DEG has given rise to a host of new 

n n 

and unusual phenomena such as the anomalous quantum Hall effect |2[ with profound 
effects on transport in these systems. The transport properties of graphene are currently 

n 

being explored in the presence of nonuniform potentials, such as in p — n junctions[3], as 
well as in periodic potentials. Effects of periodic potential on electron transport in 2D 
electron systems has been the subject of continued interest, where electrical modulation of 
the 2D system can be carried out by depositing an array of parallel metallic strips on the 
surface or through two interfering laser beams 4]. More recently in graphene, electrostatic ^ 
and magnetic |6|] periodic potentials have been shown to modulate its electronic structure in 
unique ways leading to fascinating physics and possible applications. Periodic potentials are 
induced in graphene by interaction with a substrate or controlled adatom deposition j^. 
In addition, it was shown that periodic ripples in suspended graphene also induces a pe- 
riodic potential in a perpendicular electric fieldjol. Epitaxial growth of graphene on top 
of a prepatterned substrate is also a possible route to modulation of the potential seen by 
the electrons. In this work, we complement these recent studies to discuss the effects of a 
weak electric modulation on the electrical conductivity in a graphene monolayer subjected 
to an external magnetic field perpendicular to the graphene plane. Electric modulation in- 
troduces a new length scale, period of modulation, in the system giving rise to interesting 
physical effects on the transport response. Commensurabiliy (Weiss) oscillations, in addition 
to Shubnikov de Hass (SdH) oscillations, are found to occur as a result of commensurabil- 
ity of the electron cyclotron diameter at the Fermi energy and the period of the electric 
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modulation. In [10| , on the same subject, diffusive contribution to magnetoconductivity 
was considered whereas in the present work we determine colhsional and Hall contibutions 
as well. This makes this paper a complete study of electric modulation induced effects on 
electrical conductivities/resistivities in a graphene monolayer in the presence of a magnetic 
field. 

In the next section, we present the formulation of the problem and derive expressions for 
electrical conductivities in a graphene monolayer. In section III, results of numerical work 
are presented and discussed, followed by the conclusions in section IV. 



II. FORMULATION 

We consider a graphene sheet in the x — y plane. The magnectic field B is applied along 
the z— direction. The system is also subjected to a ID weak periodic modulation U{x) in 
the X— direction. The one electron Hamiltonian reads 

H = VF(T.{p + eA) + U{x) (1) 

where p is the momentum operator, a = {o"x;Cr»/} are Pauli matrices and vf{^ lO^m/s) 
characterizes the electron velocity in graphene. In the absence of modulation, i.e. for 
U{x) = and for the vector potential chosen in the Landau guage A = {0,Bx,0), the 
normalized eigenfunctions of Eq. are given by -7==( , *'^"!a>4<i!) where 0„(a;) and 
are the harmonic oscillator wavefunctions centred at Xo = Pky. n is the Landau level index, 
/ = the magnetic length and Ly the length of 2D graphene system in the y direction. 
The corresponding eigenvalue is En = fvjjgy/n where ujg = vp \/2eB/h = vpV^/l- 

The modulation potential is approximated by the first Fourier component of the periodic 
potential U{x) = VqCosKx where K = 27t /a, a is the period of modulation and is the 
constant modulation amplitude. This potential lifts the degeneracy of Landau Levels (LLs) 
and the energy becomes dependent on the position Xo of the guiding centre. Thus energy 
eigenvalues for weak modulation {Vo ^ Ep), using first order perturbation theory, are 

En,ky = En + Vn,B COS KXo (2) 

where Vn^B = ^e~'^/'^ [Ln (u) + (u)] with L„ (u) , (u) the Laguerre polynomials 
and u = We note that the electric modulation induced broadening of the energy 
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spectrum is nonuniform. The Landau bandwidth ~ Vn,B oscillates as a function of n since 
Ln (u) are oscillatory functions of index n. Vn,B at the Fermi energy can be approximated, 
using an asymptotic expression for n ^ 1 appropriate for low magnetic-field range relevant 
to the present study, as 



Vb = V^xI^cos{KR,-^) (3) 



where Rc = kpl"^ is the classical cyclotron orbit, kp = ^I'nue and Uf. is the electron number 
density. The above expression shows that Vb oscillates with i?, through i?c, and the width 
of Landau bands 2 | | becomes maximum at 

— = r^\ (^=1,2,3,...) (4) 

and vanishes at 

^ = ^-- = 1,2,3,...). (5) 

which is termed the flat band condition. The oscillations of the Landau bandwith is the 
origin of the commensurability (Weiss) oscillations and, at the same time, are responsible for 
the modulation of the amplitude and the phase of the Shubnikov-de Hass (SdH) oscillations. 

To calculate the electrical conductivity in the presence of weak modulation we use Kubo 
formula [11]. The diffusive contribution to conductivity which arises due to the scattering 
induced migration of the Larmor circle center has already been determined for a graphene 
monolayer in [l^. Our focus, in this work, will be the calculation of the coUisional contri- 
bution to the conductivity and the the Hall conductivity. 

1. COLLISION AL CONDUCTIVITY: 

To obtain collisional contribution to conductivity, we assume that electrons are elastically 
scattered by randomly distributed charge impurities as it has been shown that charged 



impurities play a key role in the transport properties of graphene near the Dirac point 



17 



18| 



This type of scattering is dominant at low temperature. The collisional conductivity when 



spin degeneracy is considered is given by [11| 



= - - «l? (6) 
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where = [exp( '^^^ + 1)]~^ is the Fermi Dirac distribution function with = f^i for 
elastic scattering, fc^ is the Boltzmann constant and fi the chemical potential. W^^r is the 
transmission rate between the one-electron states |^) and |^'), Q the volume of the system, 
e the electron charge, t{E) the relaxation time and a| = (.^| |.^) the mean value of the x 
component of the position operator when the electron is in state |^). 

Collisional conductivity arises as a result of migration of the cyclotron orbit due to scat- 
tering by charge impurities. The scattering rate W^^r is given by 

W^^' = E l^^l' I ^"'-^^'''^ - E^')- (7) 

The Fourier transform of the screened impurity potential is Uq = lixe^ j + /c^, where 
r and R are the position of electron and of impurity respectively; is the screening wave 
vector, e is the dielectric constant of the material. By performing an average over random 
distribution of impurities, {Nj = impurity density), the contribution of the unperturbed 
part of the wavefunction,|^) = |n, ky), to the scattering rate is 

= ^Elf^'^l' K^'^^l^^'-^^"''^ \n',k'^)f6{E^,,^-E^,,,,) (8) 



with 



and 



q 



|(n, ky\ e^''-^^-^) \n', k'y)f = i [J„,„,(7) + J„„i,„,^i(7)] (9) 



l^n,n'(7)l' = ^Tj^-V-' [L:7' (7)] ; n' < n. (10) 
q\ = {ql + qy). Inserting Eq. (IHl) in Eq. ([6]) we obtain 



Here Ao = L^Ly is the area of the graphene monolayer and 7 = I {ql + Qy)/'^ = with 



Aq AqH 



(11) 

2tt 00 

with fn^ky = fiEnMy), the Fermi Dirac distribution function. Taking ^ jdipjd'y 

q 

and qy = q± simp, \Uqf ~ \Uof in Eq. f[TT]) .we obtain 



.col _ e'PN 



"^"^ A„h 



- Pof E fn,kyil - fn,ky) j^l[Jn,n'{l) + Jn~l,n'-l{l)f d-^^iE^^ky " K',fcJ- 
n,n',ky g 

(12) 



Using the following integral identity ll|, ll6| : 



7 [Jn,n'{l)] dl = j le-' [Ln{l)Y dl = {2n + 1) (13) 



where for n = n', [J„,„'(7)]^ = e~'^ [L„(7)]^ with the result 

oo oo 

j 7 [Jn-i.n'~i{l)f dj = j 7e-^ [Ln-i{l)f ^7 = (2n - 1) (14) 



oo 

lJn,n'{l)Jn-l,n'-lil)d'J = j 'jc"^ [Lnil)] [L„_i(7)] d'J = 0. (15) 



Finally, replacing the 6 function by a Lorentzian of zero shift and constant width T, ^ — »• 

ky 

J dky, y4o — ^ L^Ly, and performing the sum on n', keeping only the dominant term n' = n 
^ 

in Eq. f|T2l) . we obtain the following result 



a/l 

n=0 



OO ' „ 

dkypfn,kyil-fn,ky)- (16) 





2. DIFFUSIVE CONDUCTIVITY: 

For completeness, we also present the result for diffusive conductivity which was deter- 
mined in[10|, 

n=0 

where r is the constant scattering time and ~^q^^ = exp(3{E — i?i?)/[exp I3{E — Ep) + 1]^. 
Now a = CT^°^ + rr'^^f-f 

iNUW Uyy U^^ T Uyy . 



3. HALL CONDUCTIVITY: 

The nondiagonal contribution to conductivity [11"] is given by 

= - h') (ei vy (e'l \i) -— — ^. (18) 

{E^-E^>) 



Since /5(l-/^0(l-e 



f\i{l — f\) and VL ^ Ao = L^Ly, we obtain 



a. 



yx 



Since the x and y components of velocity operator are = and Vy 
Ho = Vpcj-ip + eA). Therefore, = vpCTx and Vy = vpcTy. Hence 



(19) 



^ when 

dpy 



(20) 



and 



(0 10) = ky\ ^ K) = '"F- (21) 
Substituting the values of the matrix elements of velocity in Eq. (|T9|) yields 



(^yx 



y 



\2 • 



(22) 



Since 



En + Vn^BCOiiKxo whcrc En = hjjgy/n and Vn^s 



^e-"/2 [Ln (u) + Ln-i (u)] we obtain 



{E^ — E^i)^ = f^^uj"] ^/n + 1 — \fn + A„ cos Kx^ 
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where 



A. 



(23) 



(24) 



Substituting Eq. (123|) in Eq. (1221) we obtain the Hall conductivity in graphene as 



a. 



yx 



h a 



a/l^ 



n=0 



J2 dk 



fn,ky fn+l,ky 



[\/n + 1 — ^Jn + A„ COS Kx^ ' 



(25) 



Elements of the resistivity tensor p^,^(/i,z/=x,?/) can be determined from those of the con- 
ductivity tensor Op^y, obtained above, using the expressions: pxx = <7yy /S, pyy = o^x / S and 
Pxy = -CTyx /S where S = a^x Oyy- a^y Oy^ with S al.y = n\e^ jB"^. 



III. RESULTS AND DISCUSSION 

The above expressions for the (coUisional, diffusive and Hall) conductivities, Eqs. (ITBI) . 
(17) and (!25|) are the principal results of this work. The integrals appearing in these equa- 
tions are evaluated numerically and the results are presented in Figure (la) at temperature 
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T = 2 K ioT a. graphene monolayer with electron density rie = 3.0 x lO^^cm ^, electric 
modulation strength Vo = 0.5meV with period a = 350nm. In addition, the following 



parameters were employed [18|, Il9|, |20|: r = 4 X 10 ^'^ s, F = OAmeV, impurity density 



Nj = 2.5 X lO^^cm^^ and e = 3.9 (using Si02 as the substrate material). We observe that 
SdH oscillations are visible in collisional conductivity a^^ whereas the Hall conductivity 
decreases with increasing magnetic field, B. Furthermore, Weiss oscillations superimposed 
on SdH oscillations are seen in ayy. To highlight the effects of modulation, we also calcu- 
late the correction to the conductivity (change in conductivity) as a result of modulation 
which is expressed as Acr^^ = cr^i^iVo) — cr^uiYo = 0) and is shown in Figure (lb). Electric 
modulation acting on the system results in a positive contribution to Acyy and a negative 
contribution to Aa^x whereas Aay^ oscillates around zero. We find that Aayy ^ Aa^x, 
which is a consequence of the fact that Aaxx has only collisional contribution, while Aayy, 
in addition to the collisional part, has contributions due to band conduction which are much 
larger. It is also seen that the oscillations in Aaxx and Aayy are 180° out of phase. To de- 
termine the effects of temperature on magnetoconductivities, comparison of conductivities 
and corrections to the conductivities at two different temperatures T = 2K (solid curve) 
and T = 6K (broken curve) are presented in Figures (2) and (3) respectively. Aaxx shows 
strong temperature dependence which is a clear signature that SdH oscillations are domi- 
nant here. Oscillations in Aayy show comparitively weaker dependence on temperature as 
Weiss oscillations, that are weakly dependent on temperature, play a more significant role 
in ayy. Furthermore, Weiss oscillations are also seen in Aayx and they are weakly sensitive 
to temperature a low magnetic fields (that is when B < 0.188T). In graphene system, the 
value of B defining the boundary between SdH and Weiss oscillations is quite low (it lies 
between 0.1 and 0.15Tesla). For smaller values of B, the amplitude of Weiss oscillations 
remain essentially the same at various temperatures. When B is large, SdH oscillations dom- 
inate and the amplitude of oscillations gets reduced considerably at comparatively higher 
temperatures. However, oscillatory phenomenon still persists. 

It can be seen from Figure (la), (2a) and (2b) that amplitude of SdH oscillations re- 
mains large at those values of the magnetic field where the fiat band condition is satis- 
fied i.e at B{Tesla) = 0.6897,0.2956,0.1881,0.1379,0.1089... when i = 1,2, 3, 4.. .in Eq.([5]) 
while is supressed at the maximum bandwidth/broad band condition, i.e at B{Tesla) = 
0.4138,0.2299,0.1592,0.1217,0.0985,... for i = 1,2,3,4,... in Eq.l^. Furthermore, zeros in 
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Acr^i, appear in close agreement with values predicted from the flat band condition. The 
amplitude of Aaxx and Acx^y becomes maximum at the broad band condition(as seen in Fig- 
ure (3)), whereas the amplitude of Acr^a; crosses the zero level at the broad band conditon 
and than a phase change of amplitude occurs. 

Components of the resistivity tensor p^^ have also been computed and shown in Figure 
(4a) as a function of B for T = 2K (solid curve) and QK (broken curve) respectively. The 
correction (change) in p^^j due to the modulation is shown in Figure (4b). To verify our 
results, we compare them in the absence of modulation with the unmodulated experimental 
results presented in [l4]. In order to carry this out, we note that the number density rie is 
related to the gate voltage (V^) through the relationship [3] = eotVg/te, where eg and e are 
the permitivities for free space and the dielectric constant of graphene, respectively, e is the 
electron charge and t the thickness of the sample. It yields Vg = A.8V for rig = 3.0 x lO^^cm"^. 
We find that the results for magnetoresistivities obtained in this work are in good agreement 
with the values given in reference [l^ for the unmodulated case at Vg = 4.8V. 

We observe in Figure (4), that the dominant effect of Weiss oscillations appears in p^x 
as it is proportional to Cyy whereas the amplitude of oscillations in pyy show a monotonic 
increase in ampitude with magnetic field signifying dominance of SdH in pyy. In Figure (5), 
we observe that the oscillations in Apxx and Apyy are out of phase and the amplitude of the 
oscillation in Apxx is greater than the amplitude of oscillation in Apyy. The out of phase 
character of the oscillations can be understood by realizing that the conduction along the 
modulation direction, which contributes to pyy, occurs due to hopping between Landau states 
and it is minimum when the density of states at the Fermi level is minimum. Oscillations in 
Pxx are much larger than those in pyy as a new mechanism of conduction due to modulation 
contributes to pxx- To highlight temperature effects on the modulated system, we present 
in Figure (6), corrections to magnetoresistivities at two different temperatures {2K, solid 
curve and QK, broken curve). These results exhibit SdH oscillation when B becomes greater 
than 0.188T as seen in Figures (5) and (6). The Weiss oscillations in Apxx are in phase 
with those of Apxy. From Figure (5c) one might infer that Hall resistivity is not affected 
by modulation. This is not so, as even Hall resistivity carries modulation effects and that is 
seen if we draw the slope of pxy as a function of magnetic field (Figure (7)). 

In order to quantatively analyze the results presented in the figures we consider the density 
of states (DOS) of this system. At finite temperature, the oscillatory part of resistivities 
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[Ap/po) are proportional to the oscillatory part of the density of states (DOS) at the Fermi 



energy, A{T /T,)AD{Ef) / where A{T/T, 



)/ sinhf;^), Do is the DOS and po is 



resistivity in the absence of magnetic field, respectively[l5']. For not too small magnetic 

2 2 

fields {B > 0.05T ), Ap/po — {ugrY Aa / ao to a good approximation, where cTq = ^-^tDo 
represents conductivity at zero magnetic field and r is the relaxation time. The analytic 
expression for the density of states (DOS) of a graphene monolayer in the presence of a 
magnetic field subjected to electric modulation has been derived in the Appendix. The 
DOS at energy E is given as 

^ 27r 

1 + 2 > — / cos[27rA;(e — t>B cost)](it exp(— 27r/c77) 
r^27r J 



D{E,Vb)=Do 



k=l 



(26) 



Do 



1 + 2^ cos(27rte) Jo(2vrA;f b) exp(— 27rA;?7) 



k=l 



where Do 



2E 



FE 



and vb = ^^Y' '^o(^) Bessel 



function of order zero. Since exp(— 27rA;?7) ^ 1 for weak magnetic fields, it is usually a good 
approximation to keep only the k = 1 term in the sum: D{E) c:^ Do + ADi{E) with 

ADi{E) 



D^ 



2 cos{2iTe)Jo{2TTVB) exp(— 2717]). 



(27) 



To determine the effects of an external magnetic field on the conductivities/resistivities of the 
system we consider Eq.(26). With a decrease in B, vb oscillates periodically with respect to 
1/B around vb = 0, increasing its amplitude proportionaly to \l\fB [Eq. (JSj)]. The function 
JoC^Trvs) decreases from 1 with an increase of | vb \= 0.3827 ~ 3/8 and than changes its 
sign. Therefore the oscillations of ADi{E) takes a minimum amplitude at the maximum 
bandwith conditions while | vb \ stays less than 3/8; it disappears when a maximum of | t;^ | 
touches at ~ 3/8; it reappears with an inverted sign for | vb \ larger than 3/8. Therefore, 
if we assume that Ap/[poA(T/Tc)] oc AD^Ep) / Do holds, we can find the position where 
oscillations of Ap/[poA{T /Tc)] vanish . That occurs at | Vb |= 0.19135(/ict;g)^/-E'i?. 

We can also find the period of oscillations in conductivities/resistivities from Eq. fl26|) 
as follows. We have D{E,Vb) ~ Do{l + 2 cos(27r£:) Jo(27rt;B) exp(-27rA;r7)} Do{l + 
2 cos(27r£:)(l — 7r^t>|) exp(— 27r/c?7)}. Since oc cos^(i^'-Rc — f )• The period of oscillation can 
be estimated by equating the increment of the cosine argument with vr. 



KA{R,) = TT, 
10 



(28) 



which leads to 



A 



- 



. , r^- (29) 

In our work ( ne = 3.0 x lO^^cm"^ and a = 350nm), therefore the period of oscillations 
comes out to be 1.933T ~^ which is in good agreement with the results shown in the figures. 

Damping of these oscillations with temperature can also be discussed. In Ref. [lO|, the 
temperature scale for damping of Weiss oscillations is given by KbT^^^^^ = thvp/Air'^a 
where b = {a/iy and vp = whence the result 



KbT, 



Weiss 



(30) 



To determine the damping temperature for SdH oscillations we, following Ref. [10| and 
use asymptotic expression for magnetoconductivity. For this, we use DOS (Eq. [26]) . 



Ill, 



D{E) 



2E 



1 + 2exp(-27rr7) cos(2 



TX 



+ 



(31) 



In the asymptotic limit of weak magnetic fields when many filled Landau levels occur, we take 



Ln-i and replace e "^^L^ by 1/ \J t{ ^/rm cos{2^/ml — tx / A) and inserting the continuum 

CXD 00 

approximation X] / dED{E) in Eq. [T71 we obtain the following result 



n=0 



a Ep{hujg) 



F + 2 exp{-2nr])A{T/Tf'^^) cos ( 27r 



El 



COS 



V2KI 



Ep 



where F 



of Weiss oscillations and A{T/T^ 



SdH\ 



(32) 

is the contribution 

[A7r^EpKBT/{hujgf]/smh[4TT'^EpKBT/{hujgf] is 



1 - AiT/T^"'"') + 2A{T/T^^''') cos^ (^Kl^ - \ 



the amplitude of the SdH oscillations. Therefore, the characteristic temperature of SdH 
oscillations is given by 



KbT, 



SdH 



{hhJg 



hoJn 



1 



(33) 



The amplitude of oscillations is given by A 



sinh(a;) 



where x = The amphtude of Weiss 



oscillations at -B = 0.3T are 0.9993 and 0.9938 at T = 2K and T = 6K, respectively The 
corresponding amplitudes for SdH oscillations are 0.6878 and 0.0882. The SdH amplitude 
decreases by ~ 87 percent whereas the amplitude of Weiss ocillations decreases by ~ 0.55 
percent for AK change in temperature. In Figures (2), (3), (5) and (6); the SdH amplitude 
decreases by ~ 77 percent when temperature is changed from T = 2K to T = 6K, and it 
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is in good agreement with the results obtained from Eqs. (130|) and (l33l) . It is due to the 
fact KbT^^''' > KbT^'^^ that the Weiss oscillations are more robust against temperatue 
changes. 

Finally, we compare the results obtained for the conductivity/ resistivity of graphene 
with those of a 2DEG given in pjj] . The characteristic damping temperatures for Weiss and 
SdH oscillations in 2DEG are KbT^deg = |^(^) and KbT^J^eg = respectively. In 
contrast, the corresponding damping temperatures in graphene are given by Eqs (30) and 
(33). On comparing the two temperature scales, we find that the damping temperatures 
of both oscillations in graphene are higher than that of a 2DEG. The ratio is found to 
be 7^ = "^.""^ ^ 4.2: where m* is the electron mass in a 2DEG and T^e is the critical 
temperature of a 2DEG; which implies that a comparatively higher temperature is required 
for damping of oscillations in graphene. This is due to the higher Fermi velocity of Dirac 
electrons in graphene compared to standard electrons in a 2DEG systems. It is evident from 
the numerical results that both, Sdh and Weiss-type oscillations, are more enhanced and 
more robust against temperature in graphene. 

To conclude, we have investigated the effects of a weak periodic electric modulation on 
the conductivity of a graphene monolayer subjected to a perpendicular magnetic field. As a 
result of modulation a new length scale, period of modulation, enters the system leading to 
commensurability oscillations in the diffusive, coUisional and Hall contributions to conduc- 
tivities/resistivities. These modulation induced effects on graphene magnetotransport are 
discussed in detail in this work. 

IV. APPENDIX 

Here we derive the expression for the density of states, Eq. (1261) in the text. We consider 
monolayer graphene subjected to a uniform quantizing magnetic field B = Bz in the presence 
of an additional weak periodic modulation potential. The energy spectrum in the quasi 
classical approximation, i.e. when many Landau bands are filled may be written as 

En,Xo = VnhWg + Vn,B COS KXo (34) 

where Vn^B = Xfe~'^/'^ [Ln (u) + Ln~i (u)] with L„ (u) , (u) the Laguerre polynomials and 
u = K'^P/2. For large n; Ln (u) ~ Ln-i (u) and Vn,B = Voe~'^^'^Ln (u). Using the asymptotic 
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expression for the Laguerre polynomials[16!]; e (u) . ^ cos{2y/nu — j) and taking 



the continuum hmit n — > |(:;^)^, where Vp = iOgl/^/2 we get 

K.^^„.-'^(iA-/^^)""co=(v2«^-^) (35) 



1/4 / 771 



/I E / 



4 



(36) 



To obtain a more general result which will lead to the result that we require as a limiting 
case we consider impurity broadened Landau levels. The self energy may be expressed as 

^'(^)-r;E/^ ,_,J_,_,,, (37, 

n Q 

which yields 

^ " / ~^E-i:~{E)- Vn,B COS Kxo - ^nwg ■ ^^^^ 



— oo 



To is the broadening of the levels due to the presence of impurities. The density of states is 
related to the self energy through 



D{E) = Im 



7r2/2r2 



(39) 



oo 



The residue theorem has been used to sum the series ^/(ri) = — {Sumof residues of 

, , — oo 

n oo ^ 

7r(cot7™)/(n)at all poles of /(n)}[l2]. Here f{n) = ^^z^ with b = Tl, c = E - 

— oo 

Ti^{E) — Vn^B COS Kxo and d = hujg. The function f{n) has a pole at c^/rf^ and the residue of 

2 oo 2 

(7r(cot 7r?2)/(n)) at the pole is ^27rcot(^). Hence ^f{n) = ^7rcot(^) and we obtain 

— oo 



y ^cot i^j^lE - 2{S-(E) + K,BCOs(Kxo)}]j 



(40) 





Separating S~(£:) into real and imaginary parts 



S-(E) = A(E) + z^^, (41) 
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Eq. (Ho!) takes the form 

^(E) f^^^) - ^""^"^ ?rfXo Sin2n + 2sinh2i; 
^ ^ 2 {hujgy J a cosh2w-cos2M ^ 



where 

M = [e - 2{ A(^) + K.B cos(i^a;o) }] (43) 

. = ^-^^ (44) 

a " oo 

27iT'^E fdxo sinh2t> 2'7iTlE fdxo,^ x / ^, x 

[^"(^'1 = (fi^^F y - co.h2„-co.2. = / -(l + 2E-(2*-)-P(-2M. 

'^-^ 

(45) 

If we define dimensionless variables e = (7^)^, V = and vb = j^yi the density of 

states is obtained as 

D{E,Vb) = Do{E){1 + 2Y^ / — cos[27Tk{e - Vb cos Kox)]exp{-27Tkr])} (46) 
where Do{E) = jj-^^YttP' = t in the above expression results in 

1 /■ 

D{E,Vb) = Do(E){1 + 2V— / cos[27rA;(e- T;BCOst)]c/texp(-27rA;?7)}. (47) 

fc=i \ 

Solving the integeral yields 

00 

D{E, Vb) = DoiE){l + 2^ cos(27rfe) Jo(27rA;t;B) exp(-27rA;?7)}. (48) 

k=l 

* Corresponding author: ksabeeh@qau.edu. pk 
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